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Abstract These notes introduce the fundamentals of black hole geometry, the ther- 
mality of the vacuum, and the Hawking effect, in spacetime and its analogues. Stim- 
ulated emission of Hawking radiation, the trans-Planckian question, short wave- 
length dispersion, and white hole radiation in the setting of analogue models are 
also discussed. No prior knowledge of differential geometry, general relativity, or 
quantum field theory in curved spacetime is assumed. The discussion attempts to 
capture the essence of these topics without oversimplification. 
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Fig. 1 The light cone at an event p. The event A is future timelike related to p, while B, C, D, and 
E respectively are future lightlike, spacelike, past lightlike, and past timelike related to p. 



1 Spacetime geometry and black holes 

In this section I explain how black holes are described in general relativity, starting 
with the example of a spherical black hole, and followed by the 1+1 dimensional 
generaUzation that figures in many analogue models. Next I discuss how symmetries 
and conservation laws are formulated in this setting, and how negative energy states 
arise. Finally, I introduce the concepts of Killing horizon and surface gravity, and 
illustrate them with the Rindler or acceleration horizon, which forms the template 
for all horizons. 



1.1 Spacetime geometry 

The line element or metric ds^ assigns a number to any infinitesimal displacement 
in spacetime. In a flat spacetime in a Minkowski coordinate system it takes the form 



where f is the time coordinate, x^y^z are the spatial Cartesian coordinates, and c is 
the speed of light. Hereafter I will mostly employ units with c—\ except when dis- 
cussing analogue models (for which c may depend on position and time when using 
the Newtonian t coordinate) . When ds^ = the displacement is called lightlike, or 
null. The set of such displacements at each event p forms a double cone with ver- 
tex at p and spherical cross sections, called the light cone or null cone (see Fig.[T]). 
Events outside the light cone are spacelike related to p, while events inside the cone 
are either /Mfwre timelike or past timelike related to p. For timelike displacements, 
ds^ determines the square of the corresponding proper time interval. 

The metric also defines the spacetime inner product g{v,w) between two 4- 
vectors v and w, that is. 



ds^ = c^dt"^ - {dx^ +dy^ +dz^), 



(1) 
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g{v,w) — ds^(v,w) — c^dt{v)dt{w) — [dx{v)dx{w) +dy{v)dy{w) +dz{v)dz{w)]. 

(2) 

Here dt{v) = v"dat = v' is the rate of change of the t coordinate along v, etc. 
In a general curved spacetime the metric takes the form 

ds^ — gfjip {x)dx"dx^ , (3) 

where {x"} are coordinates that label the points in a patch of a spacetime (perhaps 
the whole spacetime), and there is an implicit summation over the values of the in- 
dices a and /3. The metric components g^p are functions of the coordinates, denoted 
X in ([3jl. In order to define a metric with Minkowski signature, the matrix ^£,^3 must 
have one positive and three negative eigenvalues at each point. Then local inertial 
coordinates can be chosen in the neighborhood any point p such that (i) the metric 
has the Minkowski form ([TJ at /? and (ii) the first partial derivatives of the metric 
vanish at p. In two spacetime dimensions there are 9 independent second partials of 
the metric at a point. These can be modified by a change of coordinates — > jc'^ , 
but the relevant freedom resides in the third order Taylor expansion coefficients 
{d^x'^ /dx°'dx^dx'^)p, of which only 8 are independent because of the symmetry 
of mixed partials. The discrepancy 9 — 8 = 1 measures the number of independent 
second partials of the metric that cannot be set to zero at p, which is the same as 
the number of independent components of the Riemann curvature tensor at p. So a 
single curvature scalar characterizes the curvature in a two dimensional spacetime. 
In four dimensions the count is 100 — 80 = 20. 



1.2 Spherical black hole 

The Einstein equation has a unique (up to coordinate changes) spherical solution in 
vacuum for each mass, called the Schwarzschild spacetime. 

1.2.1 Schwarzschild coordinates 

The line element in so-called Schwarzschild coordinates is given by 

ds^ =[l-'j;^df -[l-'^^ ^ dr^-r^idO^ + sin^Od^^). (4) 

Here r^ = 2GM/c^ is the Schwarzschild radius, with M is the mass, and c is set to 
1 . Far from the black hole, M determines the force of attraction in the Newtonian 
limit, and Mc^ is the total energy of the spacetime. 

The spherical symmetry is manifest in the form of the line element. The coor- 
dinates 9 and (j) are standard spherical coordinates, while r measures I /In times 
the circumference of a great circle, or the square root of 1 /4-n times the area of a 
sphere. The value r = r^ corresponds to the event horizon, as will be explained, and 
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Fig. 2 Gravitational redshift. Two liglitrays propagating from r,, to rh, separated by a coordinate 
time St. Tiie corresponding proper time at r„ is less than that at . 



the value r = is the "center", where the gravitational tidal force (curvature of the 
spacetime) is infinite. Note that r should not to be thought of as the radial distance 
to r = 0. That distance isn't well defined until a spacetime path is chosen. (A path 
at constant t does not reach any r < r^.) 

The coordinate f is the Schwarzschild time. It measures proper time at r = oo, 
wheras at any other fixed r, 9, (j) the proper time interval is At = i/l — rs/rdt. 
The coefficients in the line element are independent of t, hence the spacetime has a 
symmetry under t translation. This is ordinary time translation symmetry at r — oo, 
but it becomes a lightlike translation at r = r,, and a space translation symmetry for 
r < Tj, since the coefficient of dt^ is negative there. The defining property of the 
Schwarzschild time coordinate, other than that it measures proper time in the rest 
frame of the black hole at infinity, is that surfaces of constant t are orthogonal, in the 
spacetime sense, to the direction of the time-translation symmetry, i.e. to the lines of 
constant (r, 0, 0): there are no off -diagonal terms in the line element. But this nice 
property is also why t is ill behaved at the horizon. 



Redshift and horizon 

Suppose a Ught wave is generated with coordinate period At at some radius r^, and 
propagates to another radius (see Fig. [2]l. Because of the time translation sym- 
metry of the spacetime, the coordinate period of the wave at r^, will also be At. The 
ratio of the proper time periods will thus be ATh/Ati, = ^1 — r,./r„/^l — r^./r/,, 
and the ratio of the frequencies will the the reciprocal. This is the gravitational red- 
shift. Note that as r^, the redshift is infinite. The infinite redshift surface r — r^ 
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Fig. 3 A null surface is tangent to the local light cone. 



of the spherical black hole is the (stationary) event horizon. The same is true of the 
1+1 dimensional black holes we focus on later in these notes. 

It is worth emphasizing that for a non-spherical stationary black hole, for in- 
stance a rotating black hole, the infinite redshift surface, where the time-translation 
symmetry becomes lightlike, is generally not the event horizon, because it is a time- 
like surface. A timelike surface can be crossed in either direction. In order to be a 
horizon, a surface must be tangent to the local light cone at each point, so that it 
cannot be crossed from inside to outside without going faster than light. At each 
point of such a null surface there is one null tangent direction, and all other tangent 
directions are spacelike and orthogonal to the null direction (see Fig. [3]). Therefore 
the null tangent direction is orthogonal to all directions in the surface, i.e. the null 
tangent is also the normal. If the horizon is a constant r surface, then the gradient 
V^r is also orthogonal to all directions in the surface, so it must be parallel to the 
null normal. This means that it is a null (co)vector, hence ^"'^ V^rV^gr — g" — at 
the horizon. 



1.2.2 Painleve-GuUstrand coordinates 

A new time coordinate t that is well behaved at the horizon can be defined by f = 
t+h{r) for a suitable function h{r) whose bad behavior at r^ cancels that of f. This 
property of course leaves a huge freedom in h{r), but a particularly nice choice is 
defined by 




(5) 
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where now I have adopted units with rs — 1. It is easy to see that the t-r part of the 
Schwarzschild Hne element takes the form 




ds^ =dt^-\dr+J-dt] -r^{de^- + sin^ed(j)^) (6) 



= (^l-^^dt^-^dtdr-dr^-r^{de^ + sm^ed(l)^) (7) 

The new coordinate f is called the Painleve-Gullstmnd (PG) time. At r = 1 the 
metric coefficients are all regular, and indeed the coordinates are all well behaved 
there. According to (j?]), we have ds^ = along a line of constant {r — 1,0,0), so 
such a line is lightlike. Such lines generate the event horizon of the black hole. The 
PG time coordinate has some remarkable properties: 

• the constant t surfaces are flat, Euclidean spaces; 

• the radial worldlines orthogonal to the constant t surfaces are timelike geodesies 
(free-fall trajectories) along which dt is the proper time. 

For some practice in spacetime geometry, let me take you through verifying these 
properties. Setting dt —Qin the line element we see immediately that {r, 0, 0} are 
standard spherical coordinates in Euclidean space. To find the direction orthogonal 
to a constant t surface we could note that the gradient V^f has vanishing contraction 
with any vector tangent to this surface, which implies that the contravariant vector 
g"^'^pt, formed by contraction with the inverse metric g"^, is orthogonal to the 
surface. Alternatively, we need not compute the inverse metric, since the form of 
the line element (j6]l allows us to read off the orthogonal direction "by inspection" as 
follows. Consider the inner product of two 4-vectors v and w in this metric. 



giv,w) = dt{v)dt{w) -\dr+^\ dt j (v) \^dr+^J]: dt j (w) (8) 

-p-d9 {v)dB ( w) - r2 sin^ dip {v)d(p (w) , (9) 

using the notation of Eq. (j2]). If the vector v is tangent to the constant t surface, 
then dt{v) ~ 0, so the first term vanishes. The remaining terms will vanish if 



dr + y J dt j (w) = d9{w) = d(j){w) = 0. Thus radial curves with dr- 

dO = d(j) ~ are orthogonal to the surface, and along them ds^ = dt^, i.e. dt mea- 
sures proper time along those curves. Moreover, any other time like curve connecting 
the same two spacetime points will have shorter proper time, because the negative 
terms in ds-^ will contribute. The proper time is thus stationary with respect to first 
order variations of the curve, which is the defining property of a geodesic |^ 



Even if the other terms in the line element |6|l had not been negative, they would not contribute to 
the first order variation in the proper time away from a path with (dr+^/YJrdt) =d6^ = d(p^ = 0, 
since the line element is quadratic in these terms. Thus the curve would still have been a geodesic 
(although the metric signature would not be Lorentzian). 
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Fig. 4 Painleve-Gullstrand coordinate grid for Schwarzschild black hole. Vertical lines have con- 
stant r, horizontal lines have constant t. Shown are one ingoing radial light ray and three outgoing 
ones. The one outside the horizon escapes to larger radii, the one on the horizon remains at r^, and 
the one inside the horizon falls to smaller radii and into the singularity at r = 0. 



1.2.3 Spacetime diagram of the black hole 

The nature of the unusual geometry of the black hole spacetime can be grasped 
rather easily with the aid of a spacetime diagram (see Fig.|4]l. For the Schwarzschild 
black hole, we may exploit the spherical symmetry and plot just a fixed value of 
the spherical angles (0,^), and we may plot the lines of constant r vertically and 
the lines of constant PG time t horizontally. Then the time translation symmetry 
corresponds to a vertical translation symmetry of the diagram. 

The diagram comes alive when the light cones are plotted. At a given event, the 
light cone is determined by ds^ = 0, which for radial displacements corresponds to 
the two slopes 

dt/dr= — J — (radial ligh trays) (10) 

Far from the horizon these are the outgoing and incoming lightrays dt/dr-^ ±1. 
The ingoing slope is negative and gets smaller in absolute value as r decreases, ap- 
proaching as r — > 0. The outgoing slope grows as r decreases, until reaching infin- 
ity at the horizon at r = 1. Inside the horizon it is negative, so an "outgoing" lightray 
actually propagates to smaller values of r. The outgoing slope also approaches as 
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1.2.4 Redshift of outgoing waves near the horizon 

An outgoing wave is stretched as it climbs away from the horizon. The lines of 
constant phase for an outgoing wave satisfying the relativistic wave equation are 



just the outgoing lightrays ( 10 1. The rate of change of a wavelength A is given by 
the difference of dr/dt of the lightrays on the two ends of a wavelength, hence 
dX/dt = {d/dr){dr/dt)X. The relative stretching rate is thus given by 

dX/dt d dr c 
A dr dt 2rs ' 

where in the second step the expression is evaluated at the horizon, and the dimen- 
sionful constants are restored to better illustrate the meaning. This rate is called the 
"surface gravity" K of the horizon. Later I will explain different ways in which the 
surface gravity can defined and calculated. 



We can go further and use the lightray equation ( lOi to obtain an approximate 
expression for the wave phase near the horizon. Consider an outgoing wave of the 
form e"^, with = —of + /' k{r')dr' . (This simple harmonic t dependence is exact 
because the metric is independent of t.) Along an outgoing lightray the phase is 
constant: Q — d(j) = —codt + k{r)dr, so 

^w-rAy.-^-^- (12) 

where in the second step a near horizon approximation is used, and in the last step 
the dimensionful constants are again restored. The wave thus has the near-horizon 
form 

Note that the surface gravity appears in a ratio with the wave frequency, and there is 
a logarithmic divergence in the outgoing wave phase at the horizon. 



1.3 Effective black hole and white hole spacetimes 

Many black hole analogues can be described with one spatial dimension, and I will 
focus on those here. They are simple generalizations of the radial direction for a 
spherical black hole. 

Waves or quasiparticles in a stationary 1+1 dimensional setting can often be de- 
scribed by a relativistic field in an effective spacetime defined by a metric of the 
form 

ds^ = c{xfdp- - [dx-v{x)dtf- = [c{xf - v{xf]dt^ +2v{x)dt dx- dx^ . (14) 

In fact, any stationary two dimensional metric can be put in this form, with c{x) = 1, 
by a suitable choice of coordinates (see e.g. Appendix A in Ref. f\\ for a proof 
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of this statement). If c{x) = 1 this corresponds to the PG metric, with x o r and 
v{x) ^ —l/y/r. A horizon exists in the spacetime ( 14i if |v(x)| > |c(x)| somewhere. 

The metric ( 14 1 would arise for example in a Newtonian setting of a fluid, with 
velocity v(x) in a "laboratory frame", with c{x) = c a constant speed of sound. In 
that example, the coordinate x would measure distance in the lab at fixed Newtonian 
time f, and the metric would describe the effective spacetime for waves in the fluid 
that propagate at speed c relative to the local rest frame of the fluid. If the wave 
speed in the frame of the medium depends on some ambient local conditions then 
c{x) will depend on position. 



Moving texture 

In some models the medium may be at rest in the lab, but the local conditions that 
determine the wave speed may depend on both time and space in a "texture" that 
moves. (If the motion is uniform then in the frame of the texture this is equivalent to 
the previous case.) An example of a line element of this sort is [c{y — wt)Ydt^ — dy'^. 
Here again y measures proper distance in the lab at Newtonian time f , and the texture 
moves in the y direction with constant speed w. The line element may not look 
stationary, but it has a symmetry under t t + At combined with y ^ y + wAt. 



Black hole - white hole pair 

An example that often arises has v(x) < — c(x) < in a finite interval (x_ ,x+). Then 
x+ is a black hole horizon, analogous to the one previously discussed for the PG 
spacetime, and x_ is a white hole horizon: no waves can escape from the region 
X < x+ into the region x > x+, and no waves can enter the region x > x— from the 
region x < x_. The region between the horizons is of finite size and nonsingular. 
Fig. |5]is a spacetime diagram of this scenario. Black hole horizon on the right and 
white hole horizon on the left. The vertical arrows depict the Killing vector, which 
is spacelike in the ergoregion between the horizons and timelike outside. 



1.4 Symmetries, Killing vectors, and conserved quantities 

Each symmetry of the background spacetime and fields leads to a corresponding 
conservation law. The most transparent situation is when the metric and any other 
background fields are simply independent of some coordinate. This holds for ex- 
ample with the Schwarzschild metric Q, which is independent of both t and (j) . Of 
course the spherical symmetry goes beyond just translations, but the other rota- 
tional symmetries are not manifest in this particular form of the line element. They 
could be made manifest by a change of coordinates however, but not all at once. To 
be able to talk about symmetries in a way that is independent of whether or not they 
are manifest it is useful to introduce the notion of a Killing vector field. The flow 
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of the spacetime along the integral curves of a Killing vector is a symmetry of the 
spacetime. 

Suppose translation by some particular coordinate (a indicates one particu- 
lar value of the index a) is a manifest symmetry. The metric components satisfy 
gfiv a = 0, where the comma notation denotes partial derivative with respect to x". 
The corresponding Killing vector, written in these coordinates, is — i-^- '^he 
vector with all components zero except the a component which is 1 . Then the sym- 
metry is expressed by the equation g^v.aX" = 0- This holds only in special coor- 
dinate systems adapted to the Killing vector It is not a tensor equation, since the 
partial derivative of the metric is not a tensor 

It may be helpful to understand that this condition is equivalent to the covariant, 
tensor equation for a Killing vector, 

Xw,Si+Xfi-a^Q. (15) 

where the semicolon denotes the covariant derivative. This is called Killing's equa- 
tion. One way to see the equivalence is to use the fact that in a local inertial coor- 
dinate system at a point p, the covariant derivative reduces to the partial derivative, 
and the partials of the metric are zero. Thus Killing's equation at the point p be- 
comes .^^aa +X" .ailaj} = 0, where rjoT is the Minkowski metric. This implies 
that the infinitesimal flow x*^ — > x*^ + ex"{x) generated by is, to lowest order, a 
translation plus a Lorentz transformation, i.e. a symmetry of the metric]^ 

^ For a more computational proof, note that since Killing's equation is a tensor equation it holds 
in all coordinate systems if it holds in one. In a coordinate system for which = we have 

Xa■,|}=8a^lX''■.|3=ga^l^^•|i„x'' =j{gali.a+8aa.l}-gl3cj.a)x''-ifX'' is a Killing vector the first 
tenn vanishes in this adapted coordinate system, and the remaining expression is antisymmetric in 
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A simple example is the Euclidean plane with line element ds^ = dx^ + dy^ = 
dr^ _l_ r^(i(j)2 in Cartesian and polar coordinates respectively. The rotation Killing 
vector about the origin in polar coordinates is just d^, with components 5^, as the 
metric components are independent of (j). The same Killing vector in Cartesian co- 
ordinates is xdy — ydx- This satisfies Killing's equation since Xx.x ~ = Xy.y^ ™d 

Z^,y + Zy,x = -1 + 1 =0. 



1.4.1 Ergoregions 

It is of paramount importance in black hole physics that a Killing field may be 
timelike in some regions and spacelike in other regions of a spacetime. For example 
in the Schwarzschild spacetime, say in PG coordinates (|6]l, or the 1+1 dimensional 



generalization (14 1 the Killing vector dt is timelike outside the horizon, but it is 
lightlike on the horizon and spacelike inside. For the black hole-white hole pair 
discussed above, it is the region between the black and white hole horizons (see 
Fig. [s]). This is evident because the coefficient of dt^ in the line element becomes 
negative. 

A region where an otherwise timelike Killing vector becomes spacelike is called 
an ergoregion. (The reason for the name will become clear below.) The boundary of 
this region is called the ergosurface, and it is a surface of infinite redshift, since the 
norm of the time translation Killing vector vanishes there. An ergoregion need not 
lie behind a horizon. For instance it occurs outside the horizon (as well as inside) of 
a spinning black hole. In analogue models, ergoregions can arise for example around 
a vortex |2| or in a moving soliton in superfluid ^He-A |3|. For the Schwarzschild 



black hole, and the 1+1 dimensional generalization (14i, however, the ergoregion 
always corresponds to the region inside the horizon. 



1.4.2 Conserved quantities 

Particle trajectories (both timelike and lightlike) can be determined by the varia- 
tional principle 5/ LdX = with Lagrangian L = \gpiv{x)x^x^ ■ Here A is a path 
parameter and the dot denotes dldX. The Euler-Lagrange equation is the geodesic 
equation for motion in the metric g^jy with affine parameter X. If the metric is inde- 
pendent of x" then the corresponding conjugate momentum pa — dLjdx"' — g^aX^ 
is a constant of motion. Note that this momentum can also be expressed as the inner 
product of the 4-velocity = i*' with the Killing field, u - x = gjivX^X^ ~ gfia^^ ■ 



a and j8, so adding x^;a yields zero. Conversely, if Killing's equation holds, the entire expression 
is antisymmetric in a and j8, so the first term must vanish. 
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Fig. 6 Killing energy. On the left the BCilling vector x is timelike, hence all future causal (timehke 
or lighthke) 4-momenta have positive ;(^-energy. On the right x is spacehke, hence future causal 
4-momenta like p2 can have negative ;ij-energy, while others Uke pi have positive ;|^-energy. 



Killing energy and ergoregions 

The conserved momentum conjugate to a particular timelike Killing field is called 
Killing energy. For a particle with rest mass m, the physical 4-momentum would be 
p = mu, so the Killing energy as defined above is actually the Killing energy per 
unit rest mass. For a massless particle, the physical 4-momentum is proportional to 
the lightlike 4-velocity, scaled so that the time component in a given frame is the 
energy in that frame. In both cases, the true Killing energy is the inner product of 
the 4-momentum and the Killing vector, 

^Killing ^/'■Z- (16) 

The 4-momentum of a massive particle is timelike, while that of a massless particle 
is lightlike. In both cases, for a physical state (i.e. an allowable excitation of the 
vacuum), stability of the local vacuum implies that the energy of the particle is 
positive as measured locally in any rest frame. This is equivalent to the statement 
that p is a future pointing 4-vector. 

The importance of ergoregions stems from the fact that negative Killing energy 
physical states exist there. This happens because a future pointing 4-momentum can 
of course have a negative inner product with a spacelike vector (see Fig. |6|. In an 
ergoregion, the Killing energy is what would normally be called a linear momentum 
component, and there is of course no lower limit on the linear momentum of a 
physical state. 

Penrose lIH |5l realized that the existence of an ergoregion outside a spinning 
black hole implies that energy can be extracted from the black hole by a classical 
process, at the cost of lowering the angular momentum. This is the Penrose process, 
whose existence led to the discovery of black hole thermodynamics. For a non- 
spinning black hole the ergoregion lies inside the horizon, so no classical process 
can exploit it to extract energy, but the Hawking effect is a quantum process by 
which energy is extracted. 

What do the negative Killing energy states "look Uke"? A particle with negative 
Killing energy cannot escape from the ergoregion, nor can it have fallen freely into 
the ergoregion, because Killing energy is conserved along a geodesic and it must 
have positive Killing energy if outside the ergoregion. For example, in the 1h-1 black 
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Fig. 7 Inside the horizon, the Killing vector x is spacelike, outgoing radial particles have nega- 
tive 3;;-energy, and infalling ones have positive -energy. (Since the latter come from outside the 
ergorgion, and Killing energy is conserved, they must have positive Killing energy.) 



hole, or in the radial direction of the Schwarzschild solution, a massless particle 
with negative Killing energy inside the horizon must be "outgoing" as seen by a 
local observer (see Fig.|7]i. 



1.5 Killing horizons and surface gravity 

An event horizon can be defined purely in terms of the causal structure of a space- 
time, and is meaningful even when the spacetime is not stationary, i.e. has no time 
translation symmetry. A Killing horizon on the other hand is a lightlike hypersurface 
(surface of one less dimension than the whole spacetime) generated by the flow of 
a Killing vector. This is sometimes called the horizon generating Killing vector. 

The Schwarzschild event horizon is a Killing horizon with respect to the Killing 
vector d,, as is the horizon of the 1+1 black hole. A distinction arises in the case of a 
stationary black hole with spin. Then the Killing vector dt that is a time translation 
at spatial infinity becomes lightlike at the boundary of the ergoregion, which lies 
outside the event horizon. However that boundary is timelike, so the ergosurface is 
not a Killing horizon. The event horizon of a spinning black hole is nevertheless a 
Killing horizon, but for a Killing vector di + Q^d^ that is a linear combination of 
the time translation and rotation Killing vectors, being the angular velocity of 
the horizon. In the effective spacetime of a moving texture in superfluid ^He-A, the 
horizon generating Killing vector has the similar form (9, + vv(9^, where dt and are 
time and space translation Killing vectors, and the constant w can be thought of as 
the transverse velocity of the horizon IS] . 
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Fig. 8 Boost killing flow in Minkowski space jl8| . Curves of constant (. are hyperbolic flow lines. 
Lines of constant r\ are radial from the origin, and r\ measures the hyperbolic opening angle of the 
shaded wedge. 



Rindler (acceleration) horizon 

A simple yet canonical example of a Killing horizon is the Rindler horizon in 
Minkowski spacetime. The relevant Killing symmetry here is Lorentz boosts is a 
certain direction. Geometrically, these are just hyperbolic rotations. For example, 
using the Minkowski coordinates of ([T]l a boost Killing vector is 

XB=xd,+td^. (17) 

This has covariant components {XB)a ~ "HafiXB = i^^ ^0 ™d so obviously satisfies 
Killing's equation ( [T5| . It can also be made manifest by changing from Minkowski 
to polar coordinates: 

dt^~dx^ =fdri^-df. (18) 
Then the boost symmetry is just rotation of the hyperbolic angle rj, i.e. 

Zb-^t,. (19) 

The flow lines of the Killing field are hyperbolas (see Fig. [8]l. Note that the polar 
coordinate system covers only one "Rindler wedge", e.g. jc > |f| of the Minkowski 
spacetime. The full KilUng horizon is the set |x| = |f |. 



1.5.1 Surface gravity 



Associated to a Killing horizon is a quantity K called the surface gravity. There are 
many ways to define, calculate, and think of the surface gravity. It was already intro- 
duced in Sec. 1.2.4 as the relative rate of stretching of outgoing wavelengths near 
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the horizon. I will mention here several other definitions, which are given directly 
in terms of the geometry of the horizon. 

Geometrically, the simplest definition of surface gravity may be via 

[Xla.p]X^"'% = -21^, (20) 

horizon the square bracket on indices denotes antisymmetrization, and the subscript 
H indicates that the quantity is evaluated on the horizon. That is, K is the magnitude 
of the infinitesimal Lorentz transformation generator. However the meaning of this 
is probably not very intuitive. 

The conceptually simplest definition might be the rate at which the norm of the 
Killing vector vanishes as the horizon is approached from outside. That is, 

K^Wxl.alH, (21) 

the horizon limit of the norm of the gradient of the norm of X- Notice that if the 
Killing vector is rescaled by a constant multiple x OCX^ th^n it remains a Killing 
vector, and the surface gravity for this new Killing vector is ax. This illustrates 
the important point that the intrinsic structure of a Killing horizon alone does not 
suffice to define the surface gravity. Rather, a particular normalization of the Killing 
vector is required. The symmetry implies that K is constant along a particular null 
generator of the horizon, but in general it need not be the same on all generators. 
For a discussion of conditions under which the surface gravity can be proved to be 
constant see Q. 



The surface gravity (21 1 has the interesting property that it is conformally invari- 
ant. That is, it is unchanged by a conformal rescaling of the metric gat — > ^^gah, 
provided the conformal factor Q is regular at the horizon [Tj. This follows simply 
because |x| is rescaled by £2, while the norm of its gradient is rescaled hy £2^^, and 
the contribution from d£2 vanishes since it is multipled by \x\h which vanishes. 



For the metric 



14 1 and the Killing vector x — d, we have \x\ — Vc 

1 /2 

depends on x and not t. Thus K = {—g''^dx\x\ dxlxDn ' ^'^'^ '■^^ minus sign arises 
because the gradient is spacelike outside the horizon. At a horizon where v — c 
this evaluates to \dx{v — c) \h, while at a horizon where v = — c it would instead be 
\dxiv + c)\H. 

In case c — constant, the surface gravity is thus just the gradient of the flow speed 
at the horizon. A covariant and more general version of this can be formulated. Any 
observer falling freely across a horizon can define the velocity of the static frame rel- 
ative to himself, and can evaluate the spatial gradient of this velocity in his frame. 
If he has unit Killing energy (u ■ x — 1) then it can be shown that this gradient, 
evaluated at the horizon, agrees with the surface gravity |8|. Another interesting ob- 
servation is that this velocity gradient has a sort of "cosmological" interpretation as 
the local fractional rate of expansion ("Hubble constant") of the distances separat- 
ing a family of freely falling observers stretched along the direction of the Killing 
frame velocity |8|. At the horizon, for unit energy observers, this expansion rate is 
the same as the surface gravity. 
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Computationally, a somewhat simpler definition of surface gravity is via 

[da{X^)^~2KXa]H- (22) 

This is at least well-defined: since vanishes everywhere on the Killing horizon, 
its gradient has zero contraction with all vectors tangent to the horizon. The same is 
true for Xa = gapX^^ so these two co-vectors must be parallel. If using a coordinate 
component of this equation to evaluate K, it is important that the coordinate system 



be regular at the horizon. For the metric ( 14 1, we may just evaluate the x component 
of this equation: dx{c^ — v^) = —2kXx = —^Kg^, — —2kv, which on the horizon 
V = c yields K = [dj{v — c)]h as before. (Note that this definition does not come 
with an absolute value. At a horizon v = — c it yields K = [dx{v + c)]h-) 



Surface gravity of the Rindler horizon 

The surface gravity of the Rindler horizon can be computed for example using the 



polar coordinates to evaluate (21 1. Then the norm of the Killing vector is just £, 
so (^oIxbI = which has norm 1. Thus the boost Killing vector has unit surface 
gravity. Alternatively, we may use the x component of (22 1: dxXs = -^^ — t^= 2x, 



and —2k{xb)x = 2Kt, so k = {x/t)^ = ±1 On the future horizon x — t and this is 
positive, while on the past horizon it is negative. Usually one is only interested in 
the absolute value. 

Finally, it is sometimes of interest to use the proper time along a particular hy- 
perbola rather than the hyperbolic angle as the coordinate. On the hyperbola located 
at£ — £o the proper time is c/t = iodrj . The Minkowski line element can be written 
in terms of the time coordinate T = £orj as ds^ = [l/l^Ydx^ —d£^. The scaling of 
the Killing field (9^ — {l/£o)dri that generates proper time flow on this particular 
hyperbola has surface gravity K — 1 /£q. This is also equal to the acceleration of 
the hyperbolic worldline. The relation between the surface gravity and acceleration 
can be shown quite generally using coordinate free methods, but here let's just show 
it by direct computation using Cartesian coordinates. The 4-velocity of the hyper- 
bola is the unit vector u — £Q^{x,t,0,0), and the acceleration of this worldline is 
{u-V)u^£f^^{xdt+tdx)ix,t,0,0)^£Q^{t,x,0,0).The norm of the spacelike vector 
(f ,x,0,0) is £o, so the norm of the acceleration is 1 /£o. 



2 Thermality of the vacuum 

The subject of the rest of these notes is the Hawking effect, i.e. the emission of 
thermal radiation from a black hole. The root of the Hawking effect is the thermality 
of the vacuum in flat spacetime. This thermality is known as the Unruh, or Fulling- 
Davies-Unruh, effect ||9|- In its narrowest form, this is the fact that a probe with 
uniform proper acceleration a, moving through the vacuum of a quantum field in 
flat spacetime, is thermally excited at the Unruh temperatme 



Black holes and Hawking radiation in spacetime and its analogues 



17 



Tu=ha/2nc. (23) 

(I've restored c here to show where it enters, but will immediately revert to units 
with c — \.) When described this way, however, too much attention is focused on 
the probe and its acceleration. 

Underlying the response of the probe is a rather amazing general fact: when re- 
stricted to a Rindler wedge, the vacuum of a relativistic quantum field is a canonical 
thermal state with density matrix 

p«c<exp(-2;r//^/ri), (24) 

where H-q is the "boost Hamiltonian" or "Rindler hamiltonian" generating shifts of 



the hyperbolic angle coordinate r\ defined in ( 18 1. In terms of Minkowski coordi- 
nates {t,x,y,z), Hrj is given on a f = surface of the Rindler wedge by 



H^= T,hXBdE"= xT„dxdydz, (25) 



where Tab is the energy-momentum tensor The "temperature" of the thermal state 
(|24}is 

TR=h/2n. (26) 

Like a rotation angle, the hyperbolic angle is dimensionless, so the boost generator 
and temperature have dimensions of angular momentum. 

Note that the thermal nature of the vacuum in the wedge does not refer to any par- 
ticular acceleration, and it characterizes the state even on a single time slice. Never- 
theless it does directly predict the Unruh effect. A localized probe that moves along 
a particular hyperbolic trajectory at proper distance from the vertex of the wedge 



has proper time interval dl — lodrj (cf. 18 1. When scaled to generate translations of 



this proper time the field Hamiltonian is thus ^ l^^H-q, and the corresponding 
temperature is 7b ~ Iq^Ti/Ik. The proper acceleration of that hyperbola is , so 



the probe will be excited at the Unruh temperature (23 i. 

The thermality of the vacuum in one wedge is related to entanglement between 
the quantum states in the right and left wedges. It can be understood using a simple, 
but abstract and formal, argument that employs the path integral expression for the 
ground state. Because the result is so central to the subject, I think this argument 
deserves to be explained. 

The vacuum |0) is the ground state of the field Hamiltonian H, and can therefore 
be projected out of any state \x) as |0) °^ \imt^ooe^'^\x), as long as (0|;i^) 7^ 0. The 
operator e^'^ can be thought of as the time evolution operator for an imaginary 
time —it, and its matrix elements can be represented by a path integral over fields 
on Euclidean space. This yields a path integral representation for the vacuum wave 
functional, 

fo[0]c<lim((/)|e-'^|;t)- ° ^ ^<t>e-^l\ (27) 
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Fig. 9 Euclidean Minkowski space with boundary at f = 0. When the path integral \n\ is sliced by 
constant t surfaces it presents the vacuum wave-functional. When sliced by constant angle surfaces, 
it presents matrix elements of the operator exp{ — 7tHq ), where //^ is the Lorentz boost generator. 



where S is the Euclidean action corresponding to the Hamiltonian H. The standard 
demonstration of this path integral expression for matrix elements of e^'^ proceeds 
by slicing the Euclidean space into steps of constant Euclidean time, and exploits 
the time translation invariance of the Hamiltonian. If the original Hamiltonian is 
also Lorentz boost invariant, then the Euclidean action is also rotationally invariant. 
This extra symmetry leads to an alternate interpretation of the path integral. 

Fixing a particular rotational symmetry, e.g. around the origin in the Euclidean 
tx plane, we may choose to slice the Euclidean space into steps of constant angle 
around the corresponding vertex (see Fig.|9]). This vertex divides the time slice t — 
into two halves, and the final field configuration (j) restricts to some 0^ and (j)^ on 
the left and right sides respectively. These configurations define Dirac "bras" 
and in the duals of the left and right side Hilbert spaces Jif^ and J^r. The fuU 
Hilbert space is the tensor product J^l (E) Mr. 

With this angular slicing, (and not worrying about boundary conditions at the 
vertex), we can think of the path integral as producing the matrix element of the 
operator exp(— tt//^) between 0^, regarded now as an initial state, and the final 
state 

MHM^{<^rV'''"'']\H)- (28) 

Here Hr\ is the boost Hamiltonian, which is the generator of angle shifts, and % is 
the rotation angle in the Euclidean plane. (The rotation angle is to the boost angle 
as the Euchdean time is to the Minkowski time.) The final state bra (0^1 is replaced 
by a "corresponding" initial state ket J\^l) that can be identified with a state in 
Here J — CTP^ is the operator of charge conjugation, time reversal, and reflection 
across the Rindler plane, which is a symmetry of all Lorentz invariant quantum field 
theories0 



^ For a configuration eigenstate of a real field, the ket /|(^>l) can just be identified with the same 
function (j>i, reflected by an operator P' across the Rindler plane. More generally, J includes CT to 
undo the conjugation of the (bra\ — i- \ket) duality. 
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The vacuum wave-functional ( |28] l can also be represented as a vector in the 
Hilbert space ^ J^fg, by multiplying the amplitudes ( [28] l by the corresponding 
kets and integrating over the fields: 

|0) - J^<l)L&<l>R \h)\<l>R) {<l>R\e-''"^J\<l)L) (29) 
^<j)L |fe>e-'^^')y|0^) (30) 
= ^e-'^^"|n)i|n)«. (31) 

n 

In the last line the state is expressed in terms of eigenstates \n) of the boost Hamilto- 
nian with boost energy E„ (with additional implicit quantum numbers). It is obtained 
via J\^l) = Ln-^I«)(«|0L) = L«(0z.|«)-^l«)' using the anti-linearity of 7. Then the in- 
tegral over 0i yields the identity operator, and the result follows since Hr^ commutes 
with J. The state |n) stands for the "antiparticle state" J\n). 

This exhibits the precise sense in which the quantum field degrees of freedom in 
the left and right Rindler wedges are entangled in the vacuum state. This entangle- 
ment is the origin of the correlations between the Hawking quanta and their partners, 
and it produces the entanglement entropy for quantum fields outside a horizon. Trac- 
ing over the state in the left wedge we obtain the reduced density matrix for the state 
restricted to right wedge, 

p«=Trz,|0)(0|-^e-2'^^"|«)(«|. (32) 

n 

This is the canonical thermal state (24i mentioned above The horizon entangle- 
ment entropy is the entropy of this thermal state. It diverges as the horizon area 
times the square of the momentum cutoff. 



3 Hawking effect 

The essence of the Hawking effect ifTOl is that the correlated vacuum fluctuations 
described in the previous section exist near the horizon of a black hole, which is 
locally equivalent to a Rindler horizon. The crucial difference from flat space is 
that tidal effects of curved spacetime peel apart the correlated partners. The outside 
quanta sometimes escape to infinity and sometimes fall backwards into the black 
hole, while the inside ones fall deeper into the black hole. The escaping quanta 
have a thermal spectrum with respect to the analogue of the boost Hamiltonian, 
that is, with respect to the Hamiltonian for the horizon-generating symmetry. If the 
horizon generating Killing vector is normalized to have unit surface gravity, like the 
boost KilUng vector, the temperature is again the Rindler temperature Tr ^ h/2n 



* Its matrix elements could also have been obtained directly using the wave functional |28|, via 
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( [26) l. However, for a quantum that escapes from the black hole region, the natural 
definition of energy is the generator of asymptotic time translations. For defining 
this energy we normalize the time translation Killing vector at infinity. Then the 
black hole horizon has a surface gravity K, and the temperature is the Hawking 
temperature, 

Th = hK/ln. (33) 
Note that the Unruh temperature ([23]l can be expressed in exactly the same way as 



the Hawking temperature since, as explained in Sec. 1.5.1 when the boost Killing 
field is normalized to unity on a given hyperbola the surface gravity of the Rindler 
horizon is precisely the acceleration of that hyperbola. 

For a rotating black hole, as explained in Sec. |1.5[ the horizon generating Killing 
vector is d, + Qfjd^. The eigenvalues of the Hamiltonian corresponding to this 
Killing vector ar^/s — where E and L are the energy "at infinity" and angular 
momentum respectively. Thus the Boltzmann factor for the Hawking radiation is 
^-{e-Qh^/Th ^ 'pjjg angular velocity Qh plays the role of a chemical potential for 
the angular momentum. 

Missing from this explanation of the Hawking effect is the specification of the 
incoming state. In principle, there are two places where the state can "come in" 
from: spatial infinity, and the horizon. The state coming from the horizon is deter- 
mined to be the local vacuum by a regularity condition, since anything other than 
the vacuum would be singular as a result of infinite blueshift when followed back- 
wards in time toward the horizon. This is what accounts for the universality of the 
thermal emission. However the state coming in from infinity has freedom. If it is 
the vacuum, the state is called the "Unruh state", while if it is a thermal state, as 
approprate for thermal equilibrium of a black hole with its surroundings, it is the 
"Hartle-Hawking" state. In the neighborhood of the intersection of past and future 
horizons, the Hartle-Hawking state is close to the local Minkowski vacuum. 

For black holes in general relativity, the above description of the Hawking ef- 
fect is, in a sense, the complete story. For analogue models, however, one wants a 
derivation that does not assume Lorentz invariance, and that shows the way to the 
modifications brought about by the lack thereof. Also, it is important to be able to 
allow for experimental conditions that determine different incoming states. More- 
over, in the analogue case the horizon state need not be the vacuum, since in the 
presence of Lorentz violating dispersion a different state can exist without entail- 
ing anything singular on the horizon. Thus we now take a very different viewpoint, 
analyzing the vacuum "mode by mode". It is this approach that Hawking originally 
followed when he discovered black hole radiation. It should be emphasized at the 
outset however that, unlike the previous treatment, this approach will apply only to 
free field theory, with uncoupled modes satisfying a linear field equation. 



^ The sign of the L term is opposite to that of the E term because is spaceUke while d, is timelike. 
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3.1 Mode solutions 



My aim here is to convey the essence of the Hawking effect, using a language that 
is easily adapted to analogue models in which dispersive effects play a role. Hence 
I will discuss only a system with one spatial dimension, and will highlight the role 
of the dispersion relation, using WKB methods. 

Consider a scalar field (p that satisfies the wave equation da{\J~gg"^ da(p) — 0. 



For the metric ( 14 1 we have ^J—g = c and g" = 1 /c^, g"' = v/c^, g''^ = (v^ — c^)/c^. 
Since the metric is independent of t we can find solutions with definite Killing fre- 
quency, <p = e^"'"u{x). Because of the redshift effect an outgoing solution has very 
rapid spatial oscillations of u{x) near the horizon. We can thus find an approximate 
solution near the horizon by neglecting all terms in which there is not at least one 
derivative of u{x). This yields the equation 

dx[{v^/c - c)5vm] = {2i(0v/c)d_^u. (34) 

Near a horizon x^xh where y = — c we have the expansions v / c = — I + 0{x — xh) 
and v^/c — c — —2k{x — xh) + 0[{x — xh)^] - Thus at the lowest order inx — XH the 
near horizon approximation of (|34]l becomes 



d^[{x-XH)d^u]^ {ia)/K)dxU, (35) 
whose solutions have the form 

The logarithmic divergence in the phase justifies the dominance of spatial deriva- 



tives of (p near the horizon. Note that this mode has the same form as ( 13 1, which 



we inferred in Sec. |1.2.4| using the equation of outgoing lightrays to propagate the 
phase of the wave in the near horizon region. 

Now let's see how to arrive at the same approximate solution using the dispersion 
relation with the fluid picture. First, a mode solution in a homogeneous fluid has the 
form (p ^ e^"^' e''^-" , where x is the position in the fluid frame and the dispersion 
relation is (O^ = F[k)'^ for some function F{k). For instance, for a nondispersive 
wave with speed c we have simply F{k) — ck. If the fluid is flowing with speed v 
relative to the "lab" then x = Xf + vt, where Xf is at rest with respect to the fluid. In 
terms ofxf the mode is e^''^'"^^''^^' e'^'^f , which allows us to read off the frequency as 
measured in the fluid frame, cOf ~ CO — vk. The dispersion relation holds in the fluid 
frame, so we have co — vk — ±F{k). 

If the flow velocity v{x) is not uniform, (Of = 03 — v{x)k is locally accurate pro- 
vided the change of v(x) over a wavelength is small compared to v{x) itself. The 
local dispersion relation then becomes 



(0-v{x)k^±F{k), 



(37) 
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which for a fixed Killing frequency yields a position-dependent wavevector, kio{x). 
It should be emphasized that the Killing frequency O) is a well-defined global con- 
stant for a solution, even if the Killing vector is is not everywhere timelike. 

An approximate, WKB mode solution, taking into account only the phase factor, 
is thus 

u{x)^exp^J k(o{x')dx'^ . (38) 

Finally, if the local wave velocity c{x) also depends on position in the fluid (but 
is time independent in the lab frame), then the function F{k,x) also depends on 
position. If c{x) changes slowly over a wavelength, then the mode of the same form 
is again a good approximation. For the case of relativistic dispersion F (A;, x) = c{x)k 
we obtain kg, = (o/{c + v) for the outgoing mode. Expanding around the horizon 



this yields ka){x) = {(o/k){x — xh) and so the mode takes the same form as (36 1 
derived above. 



3.2 Positive norm modes and the local vacuum 

When the field is quantized, the Hilbert space is constructed as a Fock space built 
from single particle states corresponding to (complex) solutions to the field equation 
with positive conserved "norm". The norm can be identified using a conserved inner 
product, the existence of which follows from global phase invariance of the action. 
Here I will not attempt to explain the details of this construction, which can be found 
in many expositionsPlbut instead will try to provide a simple argument that captures 
the essence of the story. In this section the relativistic case will be explained, and in 
the last section I will make some brief comments about what happens when there is 
Lorentz violating dispersion for short wavelengths. The quantum field is taken to be 
a hermitian scalar, which arises from quantization of a real scalar field. 

Positive norm modes that are localized can be recognized as those that have 
positive frequency in the fluid frame. In the relativistic case, this amounts to posi- 
tive frequency in any freely falling frame. The time derivative in the fluid frame is 



{dt)f = dt + vdx- For a mode of the form ( 36 1 near the horizon, this is dominated by 



the second term, and v sa — c, hence for such modes positive frequency with respect 
to t in the fluid frame is the same as positive frequency with respect to x. (There are 
two minus signs that cancel: v = — c < at the horizon, but the conventional def- 
inition of "positive frequency" is ^ e^"*" with o > for temporal frequency, and 
^ vvith > for spatial frequency .) 



The mode ( 36 1 with logarithmic phase divergence at the horizon can be ana- 
lytically continued across the horizon to make either a positive or a negative fre- 
quency solution. To see how this works, let's first simplify the notation a bit and 
set xh — 0, so the horizon lies at x = 0. Now a positive x-frequency function has 
the form dkf{k)e''^'^ , which is analytic in the upper-half complex x-plane since 



' For a pedagogical introduction see, e.g. 1 10|, or references therein. 
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addition of a positive imaginary part to x leaves the integral convergent. Similarly, 
a negative x-frequency function is analytic in the lower half x-plane. The argument 
of the logarithm is x = e'^ \x\, so Inx = iO + In \x\ . Continuing to —x in the upper or 
lower half plane thus gives (lnjic)± = zLin + In |x| respectively, hence 

We can thus write down positive and negative frequency continuations, 

q+=u + e-'"^/''u (40) 
q_ =e-^(o/K^^~^ (41) 

where u = 0(x)e'('"'''^''"* and u = 0(-x)e'''"/'^^'"l^l, and is a normalization factor. 
(The negative frequency continuation q_ has been multiplied by e^'^'^l'^ to better 
reflect the symmetry and thus simplify the following discussion.) 

We can now express m as a superposition of positive and negative norm parts. 

From the symmetry of the construction, the norms of q+ and q- are equal up to a 
sign, hence the ratio of the squared norms (denoted (, )) of the negative and positive 
norm parts of u is 

("+,«+) 

In the last equality I've defined the energy E = hco, and Th = hx/ln is the Hawking 
temperature. This "thermal ratio" is the signature of the Hawking effect, as indicated 
via the mode u outside the horizon. Note that this ratio is a property of the classical 
solution to the wave equation, and is determined by the ratio of the frequency to the 
surface gravity. Planck's constant enters only when we express the result in terms 
of the energy quantum hd). Note also that if the Killing vector is rescaled, then the 
Killing frequency (O and surface gravity K are rescaled in the same way, so that the 
ratio o/ff is unchanged. 

The presence of the negative frequency part m_ in u ( |42] l is unexpected from the 
WKB viewpoint. It corresponds to a negative wavevector, whereas when we solved 
the local dispersion relation we found koi){x) = {co/ k){x — xh)^^- Since the support 
of u lies outside the horizon at x > xh, it might seem that this dispersion relation 
impUes that kio (x) is positive, and thus that the frequency is purely positive. However 
this is a misconception, because a function with support on a half line cannot have 
purely positive frequency. The concept of a definite local wavevector must therefore 
have broken down. Indeed, if we examine the change of k over a wavelength we 
find {dk/dx)/k ~ (K/(o)k, which is not much smaller than k unless (O ^ K. This 
resolves the puzzleFl 

' However, it raises another one: why did the WKB type mode ~ exp{i J' kcD{x')dx') agree so well 
with the mode function l |36[ l? The answer is that \35\ is a first order equation, not a second order 
one, once an overall (5^ derivative is peeled off. 
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The local outgoing vacuum 

The local outgoing vacuum contains no outgoing excitations. More precisely, it is 
the ground state in the Fock space of outgoing positive norm modes. The outgoing 
modes we have been discussing are not themselves localized, but one can form lo- 
calized wavepackets from superpositions of them with different frequencies. Hence 
we may characterize the local outgoing vacuum by the requirement that it be anni- 
hilated by the annihilation operator^fl(^+) and a{q*_) for all positive norm modes. 
These operators can be expressed in terms of the annihilation and creation oper- 



ators corresponding to u and u using (i) linearity, (ii) equations (40i and (HTli, and 
(iii) the relation a{f) — -~a^{f*) which should be used if / has negative normrlFor 
example, a{q+) = a{u) + e^'^"'!'^ a{u) = a{u) — e^'""^'' a{u*). The vacuum condi- 
tions 

«(?+)|0)=0 (44) 
a{q*_)\0)=0 (45) 

fl(M)|0> =.e-™/''fl^(r)|0) (46) 
a(r)|0)=e-''™/''a'f(M)|0). (47) 

(48) 

If we normalize the mode u, then the commutation relation [a{u),a'' {u)] = l holds 
and implies that, in effect, a{u) — d/da^{u), and similarly for u. Thus (46i can be 



thus amount to 



solved to find the vacuum state for these particular modes of frequency O), 

|0) - exp {e--''/''a\u)aHu*)) |0z,0«), (49) 

where \QlQr) is the state with no u or ii* excitations on either side of the horizon, 
a(M)|OLOfi) = = fl(M*)|OiOfl). In flat space \QlOr) is called the (outgoing factor of 
the) "Rindler vacuum", while in a black hole spacetime it is the "Boulware vacuum". 
Expanding the exponential in (|49|) we obtain another expression for the vacuum 



|0>c<^e-"W-|„^„^)^ (50) 

where and are the number of particles in the given modej^Taking the product 
over all frequencies, we then arrive at an expression for the local vacuum of a free 
field theory near the horizon that has the same form as the general thermal result 



^ What I am calling the annihilation operator here is related to the field operator (j> by a{f) = (/, (j)), 
where / is a positive norm mode. If / is not normalized this is actually {f,f)^^^ times a true 
annihilation operator. 

' The minus sign comes from the conjugation of a factor of / in the definition of the norm, which I 
will not explain in detail here. 

Here I've use the relation (a^)"|0) = \/^.\n). 
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( [3TI obtained earlier using the path integral. The results look different only because 
here the energies of free field states with n quanta are given by ntid), and because 
here the Killing vector is not normalized to unit surface gravity. 



3.3 Stimulated emission of Hawking radiation 

So far I spoke only of the Hawking effect arising from the local vacuum at the hori- 
zon. For a real black hole this is probably the only relevant condition, but for ana- 
logue models it is possible, and even unavoidable because of thermal fluctuations, 
noise, or coherent excitations, that the in-state is not the vacuum. Then what arises 
is stimulated emission of Hawking radiation ifTTl , just as the decay of an excited 
atomic state can be stimulated by the presence of a photon. 



To quantify this process, instead of imposing the vacuum condition (44i we can 
assume the quantum field is in an excited state, 

a\q+)a{q+)\>¥) = n+\>¥) (51) 
a^(r)«(r)|f) =«-|f), (52) 

where the q± are normalized versions of ( |40|41[ ). A simple way to diagnose the 
emission is via the expectation value of the occupation number of the normalized 



mode u. Using (42 1 and ( |43] l we find 

{^'\a\u)a{u)\^') = {'¥\a\u+)a{u+)+a{u*L)a\u*L)\1') (53) 

= (M+,M+)K+e-2W^(«_ + l)] (54) 
n+ +n_ + 1 



glum I K _ 2 



(55) 



where (m_|_,m_|_) = 1/(1 — e^^Tro/ic^ Thus both n+ and n_ stimulate Hawking emis- 
sion, while only shows up in the non-thermal spectrum. Had the state been a 
coherent state, the occupation numbers would be replaced by squared amplitudes. 
Something analogueous to this occurs in the surface wave white hole radiation ex- 
periments [12] , although those waves do not have a relativistic dispersion relation. 
In the case of a Bose condensate, the appropriate in-state would presumably be more 
like a thermal state lfT3l . 



4 The trans-Planckian question 

The sonic black hole was originally conceived by Unruh fT4l in part to address what 
has come to be called the trans-Planckian question: Can the derivation of Hawking 
radiation be considered reliable given that it refers to arbitrarily high frequency field 
modes? If one assumes local Lorentz invariance at arbitrarily large boosts, then any 
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high frequency mode can be Doppler shifted to low frequency, so one might argue 
that there is nothing to be concerned about. Sometimes the point is raised that there 
is an arbitrarily large invariant center of mass energy in the collision between in- 
going and outgoing modes in the vacuum outside a horizon. However, this is true 
even in flat spacetime. We never see the effects of such collisions because they 
concern the "internal structure" of the ground state. We could presumably see this 
quantum gravity structure of the vacuum only with probes that have Planckian in- 
variant energy. Hence it is not clear to me that there is anything to worry about in 
the derivation, provided one is willing to assume local Lorentz symmetry at boost 
factors arbitrarily far beyond anything that will ever be tested. 

Even without assuming exact Lorentz symmetry, one can infer the Hawking ef- 
fect by assuming that the outgoing modes are in their local ground state near the 
horizon for free-fall frequencies high compared to, say, the light-crossing time of 
the black hole, but small compared to the Planck frequency II15I . Validity of this 
assumption is highly plausible since the black hole formation, and field propagation 
in the black hole background, is very slow compared to frequencies much higher 
than the light crossing time. One would thus expect that whatever is happening in 
the vacuum, it remains unexcited, and the outgoing modes would emerge in their 
ground state in the near horizon region. The sonic model and other analogues allow 
this hypothesis to be tested in well-understood material systems that break Lorentz 
symmetry. 

Thus one is led to consider Hawking radiation in the presence of high fre- 
quency/short wavelength dispersion, both because of the possibility that spacetime 
is Lorentz violating (LV), and because of the fact that analogue models are LV. 
However, given the very strong observational constraints on Lorentz violation lfT6ll . 
as well as the difficulty of accounting for low energy Lorentz symmetry in a the- 
ory that is LV in the UV |17|, the possibility of fundamental LV seems rather un- 
likely. Hence the main motivations for considering LV dispersion are to understand 
condensed matter analogues, and to have an example — ^probably unphysical from a 
fundamental viewpoint — in which the vacuum has strong UV modifications and the 
existence of Hawking radiation can be checked. 

The central issue in my view is the origin of the outgoing modes 1 18 1. In a con- 
densed matter model with a UV cutoff these must arise from somewhere other than 
the near horizon region, either from "superluminal" modes behind the horizon, from 
"subluminal" modes that are dragged towards the horizon and then released, or from 
no modes at all. The last scenario refers to the possibility that modes "assemble" 
from microscopic degrees of freedom in the near horizon region. This seems most 
likely the closest to what happens near a spacetime black hole, and for that reason 
deserves to be better understood. Other than a linear model that has been studied 
in the cosmological context [19|, and a linear model of quantum field theory on a 
1h-1 dimensional growing lattice ll20l . I don't know of any work focusing on how to 
characterize or study such a process. 
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5 Short wavelength dispersion 

In this concluding section, I discuss what becomes of the Hawking effect when 
the dispersion relation is Lorentz invariant ("relativistic") for long wavelengths but 
not for short wavelengths, as would be relevant for many analogue models. First I 
summarize results on the robustness of the "standard" black hole radiation spectrum, 
and then I describe the phenomena of stimulated emission and white hole radiation. 

Dispersion relations of the form — c^{k^ ±fe^/A^) have been exhaustively 
studied. The plus sign gives "superluminal" propagation at high wavevectors, while 
the minus sign gives "subluminal" propagation. Roughly speaking, a horizon (for 
long wavelengths) will emit thermal Hawking radiation in a given mode provided 
that there is a regime near the horizon in which the mode is relativistic and in the lo- 
cally defined vacuum state. This much was argued carefully in Ref. [15J, and much 
subsequent work has gone into determining the precise conditions under which this 
will happen, and the size of the deviations from the thermal spectrum, for specific 
types of dispersion relations. The dispersion determines how the outgoing modes 
arise, that is whether they come from inside or outside the horizon, and what quan- 
tum state they would be found in if the initial state were near the ground state of the 
field, as in Hawking's original calculation. 

The most recent and most complete analysis of the effects of dispersion on the 
spectrum can be found in Ref. 1211 . in which many references to earlier work can 
also be found. The basic technique used there is that of matched asymptotic ex- 
pansions, pioneered in Refs. 12211231 as applied to Hawking radiation for dispersive 
fields. The dispersive modes have associated eikonal trajectories with a turning point 
outside or inside the horizon for the sub- and super-luminal cases respectively. Away 
from the turning point approximate solutions can be found using WKB methods. If 
the background fluid velocity (or its analogue) has a linear form v(x) = —1 + Kx to 
a good approximation out beyond the turning point, then one can match a near hori- 
zon solution to WKB solutions, and use this to find the Hawking radiation state and 
correlation functions. The near horizon solution is most easily found in k space, be- 
cause while the mode equation is of higher order in x derivatives, v{x) = — 1 + iKdi^ 
is linear in k derivatives, so the mode equation is second order in d^. Further simpli- 
fications come about because a linear v{x) in fact corresponds to de Sitter spacetime, 
which has an extra symmetry that produces factorized modes. One factor is indepen- 
dent of the dispersion and has a universal (O dependence, while the other factor is 
independent of co and captures the dispersion dependence. 

The result, for dispersion relations of the form co^ — (k±k^"^^ / A^")^ (chosen 
for convenience to be a perfect square), is that the relative deviations from the ther- 
mal spectrum are no greater than of order (K'/A)(K'xiin)^('+'/^"^ times a polynomial 
in 0)/ Jfj^Here the horizon is at jc = 0, and xu,, is the largest x for which v{x) has the 
linear form to a good approximation. Thus while it is important that the Lorentz vi- 
olation wavevector scale A be much greater than the surface gravity K, this may not 



For frequencies of order the surface gravity, this quantity can also be expressed as 
(xtp/xiin)'"'"''^^", where Xtp is the (co-dependent) WKB turning point. 
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be good enough to ensure agreement with the relativistic Hawking spectrum if the 
linear regime of the velocity extends over a distance much shorter than the inverse 
surface gravity. 

At the other extreme, when the surface gravity is much larger than the largest 
frequency for which the turning point falls in the linear region, the spectrum of 
created excitations has been found to be proportional to 1 /(£>, at least for dispersion 
relations of the form O)^ — {k^ ± k'^ / A^) . This is the low frequency limit of a 
thermal spectrum, but the temperature is set not by the surface gravity but by ~ 
A (jcxiin)^/^. This result applies even in the limit of an abrupt "step" at which the 
velocity changes discontinuously from sub- to supersonic 1,24. .25 J . 



5.1 Stimulated Hawking radiation and dispersion 

For a relativistic free field, the ancestors of Hawking quanta can be traced back- 
wards in time along the horizon to the formation of the horizon, and then out to 
infinity. They are thus exponentially trans-Planckian. In the presence of dispersion, 
blueshifting is limited by the scale of dispersion, so that ancestors can be traced back 
to incoming modes with wave vectors of order A . If the dispersion is subluminal, 
those modes come from outside the black hole horizon, while if it is superluminal, 
they come from behind the horizon. Either way, they are potentially accessible to 
the control of an experiment. Instead of being in their ground state, they might be 
intentionally populated in an experiment, or they might be inadvertently thermally 
populated. Either way, they can lead to stimulated emission of Hawking radiation, 
as discussed in Section [33] 

This opportunity to probe the dependence of the emitted radiation on the incom- 
ing state is useful to experiments, and it can amplify the Hawking effect, making it 
easier to detect. Note however that when the Hawking radiation is stimulated rather 
than spontaneous, it is less quantum mechanical, and if the incoming mode is sig- 
nificantly populated it is essentially purely classical. 



5.2 White hole radiation 

A white hole is the time reverse of a black hole. Just as nothing can escape from a 
black hole horizon without going faster than light, nothing can enter a white hole 
horizon without going faster than light. Einstein's field equation is time reversal in- 
variant, so it admits white hole solutions. In fact the Schwarzschild solution is time 
reversal symmetric: when taken in its entirety it includes a white hole. A black hole 
that forms from collapse is of course not time reversal invariant, but the time re- 
verse of this spacetime is also a solution to Einstein's equation. It is not a solution 
we expect to see in Nature, however, both because we don't expect the correspond- 
ing initial condition to occur, and because, even if it did, the white hole would be 
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gravitationally unstable to forming a black hole due to accretion of matter ll26l |27]| . 
Moreover, even if there were no matter to accrete, the horizon would be classically 
and quantum mechanically unstable due to an infinite blueshift effect, as will be 
explained below. 

White hole analogues, on the other hand, can be engineered in a laboratory, and 
are amenable to experimental investigation. For example, one could be realized by 
a fluid flow with velocity decreasing from supersonic to subsonic in the direction of 
the flow. Sound waves propagating against the flow would slow down and blueshift 
as they approach the sonic point, but the blueshifting would be limited by short 
wavelength dispersion, so the white hole horizon might be stable. If the horizon is 
stable, then the time reverse of the Hawking effect will take place on a white hole 
background, and the emitted radiation will be thermal, at the Hawking temperature 
of the white hole horizon |24 | (see also Appendix D of Ref. |[T3l ). Underlying this 
relation is the fact that the modes on the white hole background are the time reverse 
of the modes on the time-reversed black hole background. Note that this means that 
the roles of the in and out modes are swapped. In particular, the incoming vacuum 
relevant to the Hawking radiation consists of low wavenumber modes propagating 
against the flow. 

When such a mode with positive norm approaches the white hole horizon, it is 
blocked and begins blueshifting. At this stage, it has become a superposition of pos- 
itive and negative co-moving frequency (and therefore positive and negative norm) 
parts. If it were relativistic at all scales, it would continue blueshifting without limit. 
It would also be unentangled with the other side of the horizon, so would evidently 
be in an excited state, not the co-moving ground state. Hence there would be a quan- 
tum instability of the vacuum in which the state becomes increasingly singular on 
the horizon. A classical perturbation would behave in a similarly unstable fashion. 

In the presence of dispersion, however, the blueshifting is arrested when the it 
reaches the dispersion scale. At that stage, if the mode becomes superluminal, it 
accelerates and both parts propagate across the horizon. If instead it becomes sub- 
luminal, then it slows down and both parts get dragged back out with the flow. In 
either case, the positive and negative norm parts are in an entangled, excited state 
that is thermal when tracing over one of the pair Thus, a dispersive wave field ex- 
hibits Hawking radiation from a white hole horizon, but with two marked differences 
when compared to black hole radiation: the Hawking quanta have high wavevectors 
even when the Hawking temperature is low, and the entangled partners propagate on 
the same side of the horizon (inside for superluminal, outside for subluminal disper- 
sion). While on the same side, the partners can separate, since in general they have 
different group velocities. 

There is an important potential complication with this story of white hole radi- 
ation. Although the singularity that would arise in the relativistic case is cured by 
dispersion, an avatar of it emerges in the form of a zero Killing frequency standing 
wave. This has been shown to arise from the zero frequency limit of the Hawking 
radiation. In that limit, the emission rate diverges as 1 /o, leading to a state with 
macroscopic occupation number that grows in time l!28ll2n . This process can also 
be seeded by classical perturbations, and it grows until nonlinear effects saturate the 
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growth. The resulting standing wave, which is a well-known phenomenon in other 
contexts, is referred to in the white hole setting as an "undulation". It is composed of 
short wavelengths that are well into the dispersive regime. Depending on the nature 
of the flow and the saturation mechanism, it could disrupt the flow and prevent a 
smooth horizon from forming. 

To conclude, I will now describe what was seen in the Vancouver experi- 
ment 1 12J. That experiment involved a flow of water in a flume tank with a velocity 
profile that produced a white hole horizon for long wavelength, shallow water, sur- 
face waves (which are dispersionless over a uniform bottom). When blueshifted 
those waves convert to deep water waves, with a lower group velocity, which be- 
have like the "subluminal" case described above. In the experiment coherent, long 
waves with nine different frequencies were launched from downstream, propagating 
back upstream towards the white hole horizon, and the resulting conversion to short 
waves was observed. The squared norm ratio of the negative and positive norm com- 
ponents of the corresponding frequency eigenmode was consistent with the thermal 
ratio (|43]){3This can be understood as coherently stimulated emission of Hawking 
radiation (see Appendix C of Ref. ITSjl for a general discussion of this process). It 
is strictly classical, but it is governed by the same mode conversion amplitudes that 
would produce spontaneous emission if the system could be prepared in the ground 
state. 



6 Acknowledgments 

I am grateful to Renaud Parentani for many helpful discussions on the material 
presented here, as well as suggestions for improving the presentation. Thanks also to 
Anton de la Fuente for helpful discussions on the path integral derivation of vacuum 
thermality. This work was supported in part by the National Science Foundation 
under Grant Nos. NSF PHY09-03572 and NSF PHYll-25915. 



References 

1. S. Corley and T. Jacobson, "Lattice black holes," Phys. Rev. D 57, 6269 (1998) [hep- 
th/9709166]. 

2. C. Barcelo, S. Liberal! and M. Visser, "analogue gravity," Living Rev. Rel. 8, 12 (2005) [Living 
Rev. Rel. 14, 3 (2011)] [gr-qc/0505065]. 



The relevant norm is determined by the action for the modes. This has been worked out assum- 
ing irrotational flow f291, which is a good approximation although the flow does develop some 
vorticity. The predicted Hawking temperature was estimated, but it is difficult to evaluate accu- 
rately because of the presence of the undulation, the fact that it depends on the flow velocity field 
that was not precisely measured, and the presence of vorticity which has not yet been included in 
an effective metric description. 



Black holes and Hawking radiation in spacetime and its analogues 



31 



3. T. A. Jacobson and G. E. Volovik, "Event horizons and ergoregions in He-3," Phys. Rev. D 58, 
064021 (1998). 

4. R. Penrose, "Gravitational collapse: The role of general relativity," Riv. Nuovo Cim. 1, 252 
(1969) [Gen. Rel. Grav. 34, 1141 (2002)]. 

5. R. Penrose and R. M. Floyd, "Extraction of rotational energy from a black hole," Nature 229, 
177 (1971). 

6. R. M. Wald, "The thermodynamics of black holes," Living Rev. Rel. 4, 6 (2001) [gr- 
qc/9912119]. 

7. T. Jacobson and G. Kang, "Conformal invariance of black hole temperatiu"e," Class. Quant. 
Grav. 10, L201 (1993) [gr-qc/9307002]. 

8. T. Jacobson and R. Parentani, "Horizon surface gravity as 2d geodesic expansion," Class. 
Quant. Grav. 25, 195009 (2008) [arXiv:0806.1677 [gr-qc]]. 

9. L. C. B. Crispino, A. Higuchi and G. E. A. Matsas, "The Unruh effect and its applications," 
Rev Mod. Phys. 80, 787 (2008) LarXiv:0710.5373 [gr-qc]]. 

10. T. Jacobson, "Introduction to quantum fields in curved space-time and the Hawking effect," 
gr-qc/0308048. 

11. R. M. Wald, "Stimulated Emission Effects in Particle Creation Near Black Holes," Phys. Rev. 
D 13, 3176(1976). 

12. S. Weinfurtner, E. W. Tedford, M. C. J. Penrice, W. G. Unruh and G. A. Lawrence, "Measure- 
ment of stimulated Hawking emission in an analogue system," Phys. Rev. Lett. 106, 021302 
(2011) [arXiv: 1008. 1911 [gr-qc]]. 

13. J. Macher and R. Parentani, "Black hole radiation in Bose-Einstein condensates," Phys. Rev. 
A 80, 043601 (2009) [arXiv:0905.3634 [cond-mat.quant-gas]]. 

14. W. G. Unruh, "Experimental black hole evaporation," Phys. Rev. Lett. 46, 1351 (1981). 

15. T. Jacobson, "Black hole radiation in the presence of a short distance cutoff," Phys. Rev. D 48, 
728 (1993) [hep-th/93031031. 

16. D. Mattingly, "Modem tests of Lorentz invariance," Living Rev. Rel. 8, 5 (2005) [gr- 
qc/0502097]. 

17. J. Collins, A. Perez and D. Sudarsky, "Lorentz invariance violation and its role in quantum 
gravity phenomenology," In *Oriti, D. (ed.): Approaches to quantum gravity* 528-547 [hep- 
th/0603002]. 

18. T. Jacobson, "On the origin of the outgoing black hole modes," Phys. Rev. D 53, 7082 (1996) 
[hep-th/9601064]. 

19. R. Parentani, "Constructing QFT's wherein Lorentz invariance is broken by dissipative effects 
in the UV," PoS QG -PR, 031 (2007) [arXiv:0709.3943 [hep-th]]. 

20. B. Z. Foster and T. Jacobson, "Quantum field theory on a growing lattice," JHEP 0408, 024 
(2004) [hep-th/0407019|. 

21. A. Coutant, R. Parentani and S. Finazzi, "Black hole radiation with short distance dispersion, 
an analytical S-matrix approach," Phys. Rev. D 85, 024021 (2012) [arXivil 108.1821 [hep-th]]. 

22. R. Brout, S. Massar, R. Parentani and P. Spindel, "Hawking radiation without transPlanckian 
frequencies," Phys. Rev. D 52, 4559 (1995) [hep-th/9506121]. 

23. S. Corley, "Computing the spectrum of black hole radiation in the presence of high frequency 
dispersion: An Analytical approach," Phys. Rev. D 57, 6280 (1998) [hep-th/9710075]. 

24. J. Macher and R. Parentani, "BlackAVhite hole radiation from dispersive theories," Phys. Rev. 
D 79, 124008 (2009) [arXiv: 0903. 2224 [hep-th]]. 

25. S. Finazzi and R. Parentani, "Hawking radiation in dispersive theories, the two regimes," Phys. 
Rev D 85, 124027 (2012) [arXiv: 1202.6015 [gr-qc]]. 

26. D. M. Eardley, "Death of White Holes in the Eariy Universe," Phys. Rev. Lett. 33, 442 (1974). 

27. C. Barrabes, R R. Brady, and E. Poisson, "Death of white holes," Phys. Rev. D 47, 2383 
(1993). 

28. C. Mayoral et al, "Acoustic white holes in flowing atomic Bose-Einstein condensates," New J. 
Phys. 13, 025007 (2011) [arXiv:1009.6196 [cond-mat.quant-gas]]. 

29. W. G. Unruh, "Irrotational, two-dimensional Siuface waves in fluids," this book, and 
arXiv:1205.6751 [gr-qc]. 



